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ARCSECOND IMAGES OF CH 3 CN TOWARD W75N 
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ABSTRACT 

CH3CN (J=6-5) was observed with a resolution of 2" toward W75N using the BIMA interferometer. 
Two continuum sources were detected at 3 mm, designated MM1 and MM2 in previous studies. Align- 
ment of two mm continuum sources with the outflow axis from MM1 suggests that these continuum 
sources may be the result of the outflow interacting with the interstellar medium. MM1 is coincident 
with compact CH3CN emission. CH3CN was not detected toward MM2. The distribution of optical 
depth (jl) is derived. An excitation analysis was not done because of large line optical depths. 



1. INTRODUCTION 

Hot cores in molecular clouds have a characteristic di- 
ameter of ~ 0.1 pc, density > 10 6 cm~ 3 and temperature 
~ 100 K. Those hot cores that are luminous in the far 
IR (Lfir > 10 4 Lq) but have no free- free radio emission 
are thought to be precursors to ultra-compact (UC) HII 
regions (Kurtz et al, 2000 and references therein). They 
are found in regions of massive star formation, often offset 
from UC HII regions (Hunter et al. 2000). Hot cores are 
principally studied through their molecular line emission 
and mm-submm continuum emission. Hofner et al. (1996) 
observed the UC HII region complex G9. 62+0. 19 and 
found that CH3CN emission is centered on the youngest, 
densest component in this complex. This component also 
drives an energetic molecular outflow (Hofner et al. 1996) 
and has weak (< 1 mJy) radio continuum emission at cm 
wavelengths (Testi et al. 2000). Based on observations of 
G9.62+0.29, Hofner et al. (1996) postulated that CH 3 CN 
emission might be a tracer of UC HII precursors. 

Wilner et al. (2001) used CH 3 CN (J=12-ll) observa- 
tions of W49N to estimate the timescale for hot core evo- 
lution; they concluded that hot cores, which preceed the 
UC HII stage, have lifetimes less than or equal to UC HII 
regions. Zhang et al. (1998) observed many molecular 
lines, including CH 3 CN (J=8-7), toward W51 and con- 
cluded from line asymmetries that CH3CN was tracing in- 
falling material. These studies found that CH3CN tends to 
be: 1) tightly confined (~ 10 4 - 10 5 AU), 2) associated with 
regions of massive star formation and 3) optically thick 
in mm-wave transitions. Watt & Mundy (1999) observed 
CH 3 CN (J=6-5) toward 4 massive star formation regions 
(2 detections, 2 nondetections). Based on emission mor- 
phology, excitation analysis and chemical models by Millar 
et al. (1997), Watt & Mundy (1999) concluded that the 
hot core near G34.26+0.15 traced by CH3CN emission was 
probably externally heated. In contrast to typical results 
from the observations cited above, they found no signifi- 



cant dust emission at 3mm. Based on single-dish detection 
and interferometer non-detection, they also concluded that 
the CH3CN emission toward Gil. 94-0. 62 has an extent of 
- 50,000 AU. 

A CH3CN (J=12-ll) survey was undertaken using the 
10m Hcinrich Hertz Submillimeter Telescope (SMT) to- 
ward 48 known massive star formation regions to de- 
termine the fraction with detectable CH3CN emission 
(Pankonin et al. 2001). Forty-six percent of the regions 
surveyed were detected in CH 3 CN emission. Because the 
angular resolution was ~30", they were unable to deter- 
mine if CH3CN originates from UC HII regions or from 
neighboring sites of even younger precursors of UC HII 
regions. Pankonin et al. (2001) found beam-averaged ro- 
tation temperatures and column densities to be consistent 
with internal heating, presumably from embedded proto- 
stars or the ionizing stars of the UC HIIs. They argued 
that beam dilution allowed only lower limits for CH3CN 
column densities. 

An investigation in the lines of CH3CN (J=6-5) using 
interferometric resolution toward W75N is presented here 
to better determine the properties of this source and to 
establish how important resolution effects may be in the 
analysis of CH3CN line emission. W75N was selected be- 
cause it is a relatively nearby region of massive star for- 
mation region with strong CH3CN emission and there are 
few previous observations of the hot core in this object. 

2. OBSERVATIONS 

Observations of W75N were obtained with the BIMA 
millimeter interferometer (Welch et al. 1996) on October 
8, 11 and 12, 1999 and April 30, 2000 in the C config- 
uration and on October 14, 1999 in the B configuration. 
Baselines ranged from 180m to 60m. The phase calibra- 
tor, 2025+337 (J2000), was observed approximately ev- 
ery 25 minutes. Our primary flux density calibrator was 
Mars, whose integrated flux density was 40 Jy on Oct. 10, 
1999. Flux calibration is accurate to 20%. 2025+337 was 
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found to have a flux density of 2.29 Jy, which is within 
the range of previous measurements (Welch et al. 1996). 
All data editing and reduction were performed using the 
MIRIAD package (Sault, Tcubcn, & Wright 1995). Nar- 
row channel observations were obtained covering the fre- 
quency range 110.3485 GHz to 110.3916 GHz covering the 
CH 3 CN J=6-5 K=0-4 emission lines. These windows were 
not wide enough to observe the nearby CH^CN transi- 
tions. Continuum subtraction was performed using two 
100 MHz windows. These two windows were inspected 
and found to be free of lines. All maps were made using 
natural weighting. The correlator setup, rms noise and 
beam sizes are summarized in Tabic 1. 

EDITOR: PLACE TABLE 1 HERE. 

3. RESULTS: CONTINUUM AND CH3CN LINE EMISSION 
3.1. MM1 

Our 3mm continuum image is shown in Figure 1 in both 
contours and greyscale. The line emission integrated over 
the K=0-4 components is shown in Figure 2 as contours 
superimposed on the 3mm continuum shown in greyscale. 
The CH 3 CN (J=6-5) spectrum integrated over a 4" x 
4" area centered on MM1 is shown in Figure 3. The 
line emission is confined to MM1. As traced by CH3CN, 
MM1 is one of the most compact hot cores yet observed 
(^deconvolved = 5.0" x 1.5" = 10,000 AU x 3,000 AU). The 
distribution of line emission is slightly elliptical, with a 
deconvolved major-to-minor axis ratio of 3.4 and position 
angle of 24° ± 10° west of north, determined by a least- 
squares 2D Gaussian fit to the line emission. The indi- 
vidual K-components were imaged separately and found 
to be indistinguishable from the integrated line emission. 
The MM1 continuum emission is also elongated with a po- 
sition angle, 24° west of north but deconvolved axis ratio 
of 1.6. MM1 has a peak 3mm continuum flux density of 80 
± 16 mJy/beam and integrated continuum flux density of 
390 ± 80 mJy (integrated over 200 square arcsecs). The 
results of our analysis are given in Table 2. 

EDITOR: PLACE FIGURE 1 HERE. 



EDITOR: PLACE FIGURE 2 HERE. 



EDITOR: PLACE FIGURE 3 HERE. 

Previous observations with the VLA show that 3 com- 
pact cm continuum sources lie within the central region of 
MM1 (Hunter et al. 1994 at 1.3 cm and 3.6cm and Tor- 
relles et al. 1997 at 1.3 cm). Those observations had beam 
sizes ~ 0.1", which permitted them to resolve the cm- 
sources (see Figure 1, each VLA source is represented by 
a white X). The sources lie roughly along the mm contin- 
uum major axis of MM1. The mm continuum and CH3CN 
major axes are both perpendicular to the MM1 outflow 
axis (discussed below). This alignment is suggestive of 
an oblate cloud core, although the mm and CH 3 CN mor- 
phology could also result from the alignment of the three 



compact VLA continuum sources. MM1 and MM2 were 
not detected by the MSX satellite (Mill et al. 1994) be- 
tween 8//m and 21/xm. Modeling the emission as thermal 
dust emission at T=35 K (see section 3.3.2), we estimate 
that expected flux densities (S„ < 1 mJy) are below the 
MSX detection limits (S„ > 100 mJy). 

3.2. MM2 

A second mm-continuum source, MM2, was found « 
5" to the southwest of MM1. MM2 has a peak 3mm 
continuum flux density of 45 ± 9 mJy/beam and inte- 
grated continuum flux density of 200 ± 40 mJy. MM2 
was not detected at cm wavelengths, indicating that its 
mm emission is probably due to dust. Previous submm 
observations do not have the spatial resolution to distin- 
guish between MM1 and MM2 (Hunter et al. 2000). MM2 
was also not detected via K-band photometry by Moore 
et al. (1988, 1991) or MSX (8/im to 21 fim). Absence of 
emission at K-band is probably due to intrinsically weak 
emission, although high optical depth may also be impor- 
tant at A < 2/im. There are four possible explanations 
for these observations which we cannot distinguish be- 
tween: 1) MM2 is heated by a proto O-star which has 
not formed a detectable HII region because of infall, 2) 
MM2 is heated by a non-ionizing star, 3) MM2 is heated 
by internal shocks, probably associated with the outflow 
from MM1, or 4) MM2 is heated externally, probably by 
the three cm-components of MM1. 

3.3. Outflow 

A bipolar outflow has been detected toward MM1 in 
CO (J=3-2) (Davis et al. 1998a; Hunter et al. 1994), CO 
(J=2-l) (Davis et al. 1998b) and in cm-radio continuum 
(Torrefies et al. 1997); all of whom find that the outflow 
position angle is 66° east of north (indicated by a black 
line in Figures 1 and 2). The outflow source has not 
been clearly determined. The outflow axis is perpendicu- 
lar to the mm-continuum and CH3CN major axis. The cm 
continuum is likely produced by free-free emission, which, 
assuming the source is optically thin, would not be de- 
tectable at 3 mm. H 2 and OH masers have also been 
measured along the outflow axis, although no significant 
velocity gradient was measured (Hunter et al. 1994 and 
Torrefies et al. 1997). Shepherd (2001) observed W75N 
at 3 mm and 1 mm and found a third mm source (MM3) 
^10" to the northeast of MM1, hinted at by the lowest 
contour in Figure 1. MM3, like MM2, lies along the out- 
flow axis but on the opposite side of MM1 relative to MM2. 
Using flux density mcasurmcnts at 1 mm (230 GHz) and 
3 mm (90 GHz) of Shepherd (2001), the spectral indices of 
MM2 and MM3 are 3.1 and 1.8, respectively. Comparing 
Shepherd's measurements with those in this study, we es- 
timate the error in the spectral index of MM2 to be large, 
-1. The alignment of MM2 and MM3 suggests that they 
are associated with the outflow from MM1. In this sce- 
nario and based on its spectral index, emission from MM2 
and MM3 could originate either from shock heated dust 
or shock-ionized optically thick gas in the outflow. 

EDITOR: PLACE TABLE 2 HERE. 
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3.4. Analysis 
3.4.1. Boltzmann Analysis and Optical Depth 

CH3CN is a symmetric top molecule whose K- 
components are closely spaced in frequency but have a 
wide range of excitation energies above ground. J=6- 
5, K=0 has an excitation energy of 13 K above ground, 
whereas the K=3 transition, which is only 19 MHz lower 
in frequency, is 76 K above ground. Assuming the line 
emission from MM1 is optically thin, an excitation anal- 
ysis, based on the radiation transfer equation, allows a 
determination of the rotation temperature (T rot ) and col- 
umn density (N tot ). 

The pattern of K-componcnts were simultaneously fit 
with Gaussian profiles. That is, we fixed the line sepa- 
rations and shifted the overall pattern to obtain a least 
squares fit. We also forced each line to have the same 
FWHM, true if all emission lines originate in the same 
volume. Thus, we allowed only 6 free parameters at each 
position: one common width (FWHM), one pattern veloc- 
ity and 4 intensities (T^). The Gaussians were fit using 
a Levenberg-Marquardt least-squares minimization tech- 
nique (Argonne National Labs, Minipak Project 1980). 
We discarded fits to spectra at any position for which 
T^ < 3er for any single emission line. 

Some of the assumptions on which the excitation analy- 
sis rests are not strictly true. Specifically, the ratios of K- 
component intensities show that tl > 1 in MM1 for many 
lines. Figure 4 shows the K=3 optical depth as calculated 
from its ratio with the K=2 line. The K=3 optical depth 
was found to range from 0.7 to ~ 8.5. If the emission lines 
originate from different volumes of gas, however, this opti- 
cal depth measurement would not be accurate. The optical 
depth of the K=0 and K=l lines were not calculated be- 
cause of blending. The optical depth of the K=4 line was 
not calculated because the line was not strong enough for 
a spatially resolved analysis. Goldsmith & Landger (1999) 
studied the effects of optical depth on the Boltzmann anal- 
ysis and concluded that in such an analysis optically thick 
emission lines can produce a large scatter. In order to 
avoid calculating erroneous excitation temperatures and 
column densities, Goldsmith & Landger (1999) recom- 
mend combining observations of lines with a wide range of 
rotational states. Pankonin et al. (2001) compared results 
between their Boltzmann analysis and three different sta- 
tistical equilibrium analyses and found that the average 
difference between the two methods in T rot was —30%. 
Pankonin et al. (2001) also derived lower limits for N tot 
which were a factor of —3 lower than the statistical equi- 
librium analyses. Because the optical depth of these lines 
is large, we limit our analysis to the determination of line 
optical depths from line data and mass from continuum 
data. Pankonin et al. (2001) found that in a sample of 20 
sources with 5 detected CH3CN (J=12-ll) K-components, 
8 had K=2 component optical depths > 8. Since W75N 
was among the 12 sources they reported without an optical 
depth measurement, our analysis here implies that CH3CN 
line emission from hot cores frequently may be optically 
thick. In such cases, a source model including kinematics 
coupled with a radiative transfer/statistical equilibrium 
analysis is essential to correctly determine physical prop- 
erties such as T r ot and N tot . Observations with signifi- 



cantly higher resolution observations, high signal-to-noise 
and a wide range of excitation energy will be required to 
provide enough constraints to justify such a model. 

The CH 3 CN emission observations toward W75N show 
some similarity to observations toward IRAS 20126+4104 
(Cesaroni et al. 1999). That study concluded that CH 3 CN 
was tracing an infalling accretion disk. The resemblance, 
however, is only strong enough to encourage further ob- 
servations as described above. 

EDITOR: PLACE FIGURE 4 HERE. 



3.4.2. Mass Estimate 

We can also approximate the mass in MM1 and MM2 
using the measured mm-continuum. We estimate the free- 
free emission contribution to the 3 mm emission by ex- 
trapolating from observations at 1.3 cm (Torrefies et al. 
1997) and 3.6 cm (Hunter et al. 1994). The three cm- 
components of MM1 all have rising spectra from 3.6 cm to 
1.3 cm, indicating the sources are partly optically thick. 
Using the spectral index between 3.6 cm and 1.3 cm to 
esimate the flux density at 3 mm due to free-free emission 
would predict ~ 100 mJy. If the sources are optically thin 
shortward of 1.3 cm, however, we estimate the free- free 
component would be ~20 mJy at A = 3 mm. Since we 
wish to obtain an upper mass limit (see below) and would 
prefer to overestimate the dust emission, we assume the 
free-free component of the 3 mm emission of MM1 is —20 
mJy. Assuming thermal dust emission accounts for the 
remaining 3 mm emission (370 mJy) and following the 
method of Hildcbrand (1983), mm-continuum emission is 
related to the total mass (gas + dust) according to: 

F„ D 2 



M, 



total 



> 



(1) 



where D is the distance to the source, F^ is the continuum 
flux density, B„ is the Planck function at frequency v and 
k v is the opacity of dust at v. M tota / is given as a lower 
limit because the measured flux density may not represent 
the entire range of dust temperatures present. We assume 

( V 

a gas-to-dust ratio of 100 and thus k v = 0.006 1 2 45GHz ) 

cm 2 g" 1 (Carpenter 2000). We assume a frequency index, 
/3 = 1.5 (Shepherd 2000, Pollack et al. 1994) for the dust 
opacity. This result is sensitive to the assumed value of 
T<j M st- If we use a lower limit for T dust = 10 K, however, 
the relationship above gives an upper mass limit of 440 
M Q for MM1 (see Table 2). If more than -500 M were 
present at any reasonable temperature, we would measure 
a greater flux density at 3 mm. Since typical values for 
the mass of clouds that give rise to O-stars are > 1,000 
Mq (see Hunter et al. 2000), we conclude that W75N is 
unlikely to be forming an O-star. This result is consis- 
tent with Hunter et al. (1994) who inferred the spectral 
types of the three cm sources comprising MM1 to be early 
B-stars. Since Hunter et al. (1994) did not account for 
dust absorption or free-free self-absorption, however, the 
true spectral types may be earlier. The mass upper limit 
calculated for MM2 using the same method is reported in 
Table 2. 
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4. CONCLUSIONS 

CH 3 CN (J=6-5) was observed toward W75N using the 
BIMA interferometer with 2" resolution. Five compact (~ 
10,000 AU) K-components were detected (K=0-4) toward 
the peak of MM1 . CH 3 CN emission is elongated with 
a deconvolved major-minor axis ratio of 3.4. The major 
axis coincides with the alignment of three cm-continuum 
sources reported previously with 0.1" resolution VLA ob- 
servations. Thus, the elongation is probably caused by 
multiple sources lying approximately along a straight line 
projected on the plane of the sky, although we cannot rule 
out CH3CN tracing out an oblate cloud core. A bipolar 
outflow was previously detected toward MM1 in CO(J=l- 
0, J=3-2) and cm-continuum. MM2 and MM3 (detected 
by Shepherd 2001) lie on opposite sides of MM1 along the 
outflow axis, indicating that they may be heated by shocks 
as the outflow interacts with the local interstellar medium. 



Analysis of K=2 and 3 emission lines indicates that 
CH 3 CN (J=6-5) is optically thick (r max w 8). Thus, a 
radiative transfer/statistical equilibrium model which en- 
corporates source kinematics and temperature and density 
gradients will be required to estimate T rot and Nch 3 cn- 
The results of this study support the previously observed 
pattern that CH3CN in massive star formation regions 
tends to be compact (~ 10 4 -10 5 AU) and optically thick 
in mm- wave transitions. 
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the Wisconsin Space Grant Consortium. E.B.C. acknowl- 
edges partial support from NSF grant AST-96 17686. We 
would like to thank an anonymous referee whose comments 
helped strengthen the paper. The discussion and conclu- 
sions in this paper do not represent the views of the NSF. 
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Table 1 
Observational Summary 



Source 
Phase center 
RA (J2000) 
DEC 

Transitions 
Center rest frequency 
Spectral resolution 
Continuum Bandwidth 
Antenna HPBW 
Synthesized HPBW 
Synthesized beam position angle 
^brightness conversion 
Max spatial scale 

HMSchannel 

RMS 

cont 

Flux calibrator 
Mars (Oct 10 1999): 
Phase calibrator 
2025+337 (Oct 8-12, 1999) 



W75N 

20 h 38 m 36.6 s 
42° 37' 32" 
CH 3 CN (J=6-5) K=0-4 
110.3701 GHz 
1.1 km/s 
200 MHz 
96" 

2.5" x 2.2" 
21° E of N 
18.2 K/ (Jy/beam) 
60" 

136 mJy/beam 
4.2 mJy/beam 

40±8 Jy 

2.3±0.4 Jy 
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Fig. 1. — 3mm continuum in grcyscale and contours. Contours arc at 2, 4, 6, 10 & 14 a (4.2 mjy/beam). H2O masers are represented by 
white plus symbols (Torrelles et al. 1997). VLA cm-radio continuum sources are represented by white Xs (Hunter ct al. 1994; Torrelles et 
al. 1997). The CO and cm-radio continuum outflow orientation is indicated by a black line (Davis et al. 1998a; Torrelles et al. 1997). The 
synthesized HPBW is indicated in the top-left corner. The origin of the outflow is not well determined. 
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Fig. 2. — CH3CN emission lines in contour over 110 GHz continuum in greyscalc. Contours arc 5, 9, 13, 17 & 21 a (8.6 mjy/bcam). The 
CO and cm-radio continuum outflow orientation is indicated by a black line (Davis ct al. 1998a; Torrelles et al. 1997). The synthesized 
HPBW is indicated in the top left corner. 
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Fig. 3. — Spectrum of CH3CN emission integrated over a 4" by 4" box centered on MM1. LSR velocity is given relative to the J=6-5 K=3 
component. The conversion to brightness temperature is 18.2 K/(Jy/beam). No reasonably abundant molecular species with low excitation 
transitions were found blended with the CH3CN emission lines. 



Table 2 

Continuum Results and Analysis 



Source 
RA 
DEC 

Peak S„ (mjy/beam) 
Integrated S„ (mJy) 
Mass upper limit (M©) 



MM1 

20 ?l 38 m 36.42 ;s ± .03 s 
42°37'34".0±0".5 
80 ± 16 
390 ± 80 
< 440 



MM2 

20 ,l 38 m 36.0 s ± .I s 
42°37'32" ± 1" 
45 ± 9 
200 ± 40 
< 240 
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Fig. 4. — Beam averaged optical depth of the K=3 component in contours overlaid upon the continuum in greyscale. Contours are 3, 5 and 
7. Maximum and minimum are 8.5 and 0.7, respectively. The optical depth was calculated from the peak T * A ratio for K=2 and K=3 lines. 



